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I. INTRODUCTION 



The QCD factorization has been widely used to investigate the charmless hadronic 

B decays. For an operator Oi of the weak effective Hamiltonian, the matrix element for 
B — > M1M2 decays under the QCD factorization is found to be expressible as 

fi 



{M,M2\0i\B) = 



duTlAu)^M2iu) + (Ml ^ M2 



+ 



(1) 



where T^'^^^^ denote the parton amplitudes and $b, $Mi, represent the light-cone dis- 
tribution amplitudes (LCDAs) for the initial state B meson and the final state Mi and M2 
mesons, respectively. The parton amplitudes T^^^^^ contain short distance interactions in- 
volved in the decay processes. The LCDAs $B) ^Nh^ and are introduced to account for 
the long distance interactions. The F^^^^{m\) with j = +, are the B Mi transition 
form factors. The meson state vector |Mj), z = 1, 2, for the meson Mi is composed of Fock 
states with different number of partons 



\Mi) = \qq)M, + \qq9)hh + 



(2) 



So far, most applications of the factorization formula Eq. ([T]) are limited to leading Fock 
state \qq)Mi of the light mesons. However, the three parton Fock state \qqg)M^ of the Mj 
meson can also contribute. 

The corrections related to the higher Fock state are usually classified as subleading twist 
contributions, since their contributions are suppressed by factors of 0{l/m^) with n > 1 in 
comparison with the leading ones. Here, mi, denotes the h quark mass. Within QCD fac- 
torization, one can employ the Feynman-diagram approach or the effective-theory approach 
for studies of subleading twist contributions. There exist established Feynman-diagram ap- 
proaches for processes other than hadronic decays, such as the calculation scheme for the 
inclusive hard scattering pro cesses or the method for the exclusive hard scatter- 



mg processes 



29 
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331 ]. However, a systematic Feynman-diagram approach for 
charmless hadronic B decays is still inaccessible. On the other hand, the effective-theory 
approaches for charmless hadronic B decays have been extensively investigated in recent 



years [34 



35|, 



36, 
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39|. 
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In this paper, a calculation scheme based on the Feynman-diagram approach for charm- 
less hadronic B decays will be developed. We will concentrate on the construction of this 
calculation scheme and apply the constructed method to calculate the tree level three par- 
ton corrections. The 0{as) three parton corrections is also desirable to understand their 
factorization properties. Since the related analysis is tedious, we plan to present the rele- 
vant calculations in our another preparing work [4^. The organization of this paper is as 
following. In Section II, the calculation scheme will be constructed. The factorization of the 
tree level three parton corrections into the partonic and hadronic parts will be outlined. In 
Section III, the analysis on how the tree level three parton corrections can be factorized into 
its partonic and hadronic parts will be described in detail. In Section IV, we will apply the 
results of the Section III to make predictions for the branching ratios of i? — > nK decays. 
The last section devotes for discussion and conclusion. 



II. COLLINEAR EXPANSION AT TREE LEVEL 

In this section, we will generalize the collinear expansion method to calculate the 

three parton corrections from the Fock state \qqg) of the meson M2 in the decay B M1M2. 
There exist other types of power corrections, such as the power corrections from soft gluons 
or renormalons. We identify these as non-partonic power corrections. For these power 
corrections, our proposed scheme may not be useful. However, to include these non-partonic 
power corrections requires further assumptions beyond the factorization. For example, the 
soft gluon power corrections are better determined by nonperturbative theories, such as the 
QCD sum rules or lattice QCD. In this work, we only investigate how the partonic (or the 
dynamic) power corrections can be included into the QCD factorization in a consistent way. 

The collinear expansion method arises from a motivation of generalizing the leading twist 
factorization theorem for the hard scattering processes to include the corrections from high 
Fock states of the target hadrons. The original idea of the collinear expansion method was 
proposed by Polizer at 1980 [4]. The systematical method was developed by Ellis, Furman- 
ski and Petronzio (EFP) in their pioneer works Using the collinear expansion, the 

EFP group showed that, for the DIS processes, the twist-4 power suppressed corrections 
can be factorized into short distance and long distance parts, which are in a similar factor- 
ized form as the leading twist contributions. However, in the EFP's approach, the parton 
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interpretation for the twist-4 corrections are lost. 







b recover the parton model picture, 



Qiu then introduced a Feynman- diagram approach [7|] to re-formula the EFP's method. In 
this Feynman- diagram language, a parton model interpretation for the twist-4 corrections 
becomes trivial. 

To begin with, we first express the matrix element of an operator Oj of the weak effective 
Hamiltonian -ffe// of the standard model for the hadronic decays B — > M1M2 in terms of 
parton model amplitudes 

{M,M2)mB) 

= E Pr'^'i^M.) [ 7^Tr[7;^(0$M.(0] + (Mi-M2) 

~^ J (27r^ {2^}^ {2n)'^ '^^ [Fi'{lN , I Ah j ^Mz ) ^b{Ib)^Mi {hh ) ^^h {h-h)] , (3) 

where p^~*^^^ denote the form factors for B —>■ Milv transition. The form factor Fj are 
defined as 

{Mr{p)\q^,h\B{P,)) = F^{q'){p + P,), + ^^^^^^^^^q, (4) 

where q = Pb — p and F^{q^) = F^i^q^) under the limit q^ 0. The parton amplitudes 
T^^j{l), Tj^-_^(/i, Z2) and T/^(/b, /mi, ^Mz) are defined to describe the hard scattering center 
involving four parton, five parton, and six parton interactions corresponding to those dia- 
grams depicted in Fig. [11(a)- (c), respectively. Note that there are also other types of parton 
amplitudes involving five or six parton interactions not being presented in Fig. [H which can 
be attributed to either the physical form factors, or to be of higher twist than three. We 
have neglected these contributions in Eq. ([3]). The Tr symbol denotes the trace operation 
applied on the color and spin indices. For convenience, we employ the light-cone coordinate 
system such that = {p'^ + q^) with two light-like vectors q'^ = (g^, q^, q]_) = {Q, 0, 0) and 
pfj. — = (0,Q,0) with Q = m^/v^, which are defined as the momenta carried 

by the final state M2 and Mi mesons, respectively. The Mi meson is defined to receive 
the spectator quark of the bottom meson. The M2 meson is defined as the emitted meson 
produced from the hard scattering center. The hadron amplitudes and $^j^(/i,Z2) are 
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defined as 

$M.(0 = J rfV(M2|g(i/)g(0)|0), (5) 

In our language, Eq. ([3]) contain leading, sub-leading and higher twist contributions. The 
Eq. ([3]) becomes meaningful only if the leading twist contributions can be separated from the 
sub-leading twist contributions. For this purpose, we employ the Qiu's Feynman- diagram 
coUinear expansion approach j3] to expand each Feynman diagram in a twist by twist man- 
ner. As the loop corrections are considered, the twist expansion then interplays with the 
expansion in a<j. To be specific, we choose the following expansion strategy. 

1. All possible Feynman diagrams ordered in are first drawn. 

2. According to the collinear expansion (developed below), each Feynman diagram of 
order 0(0;^) with n > is expanded into a series ordered by increasing twist. 

3. The contributions of the same twist order from the expansion series of each Feynman 
diagram with the same as order are added up together. 

4. The factorization properties of the final expression with a specific twist and a specific 
as order are analyzed. 

The last one is important for us to derive a meaningful perturbation theory beyond the 
leading twist. 

For latter uses, we define the soft, collinear and hard loop parton momenta. We let the 
soft momentum scale as (/"'', /~, ~ (A, A, A), the collinear momentum scale as (/"*", /~, l±) ~ 
(Q,A^/Q, A), and the hard momentum scale as ~ {Q,Q,Q)- The scale variables 

are defined as Q ~ frih and A ~ ^qcd- For a collinear loop parton, it is convenient to 
parametrize its momentum F into its components proportional to the meson momentum q^, 
the light-cone vector n'^, and the transversal directions 

F = n- + -i^jn^ + 11 , (7) 

where the vector satisfies n ■ q = 1, n ■ lj_ = 0, and = 0. For convenience, we further 
define the collinear component, F, the on-shell component, and the off-shell component, 
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= r + ^r^n'^ + , 



of the momentum as 

2n-l 

1% = -^n^ . (8) 
^ 2n-l ^ ^ 

In the above expansion of the parton momentum into different parts, we have assumed 
'^Mi = 0, i = 1,2, and = 0. The loop partons except of the bottom quark are assumed 
massless for simphcity. The contributions from non- vanishing hght quark masses are taken as 
corrections. Because the hght quark mass contributions are relatively negligible as compared 
to the bottom quark mass, we also neglect the light quark mass effects in the following 
calculations. 

According to the parametrization in Eq. ([7]), a parton propagator can be separated into 
.t. long d.tance part and .hort d.tance part (the =pecM p.opagato.). If we wnte the loop 
parton propagator as [3J 

dH 

= FL{y,z)+Fs{y,z) , (9) 

where 

.,,0.|...,0^JL. (10) 

The Fl{1) propagator corresponds to the long distance part of the propagator, since 
FL{y,z) oc Oiy — z). The Fsil) propagator represents the short distance part because 
Fsiy, z) (X 6{y — z). We now describe one important property of the long distance propaga- 
tor Fl{1). In a parton amplitude, the Fl{1) may contact with a ^n^ component of a vertex 
j^. Whenever this happens, the ^n^ vertex will extract one short distance propagator Fs{l) 
and one interaction vertex i'ji, from the relevant hadron amplitude Q] 

The momentum factor (/ — ly is then absorbed by the hadron amplitude due to the Ward 
identity We now explain how the identity Eq. (ITT!) can be obtained by a simple manipu- 
lation. We first insert an identity 1 = into the left hand side of Eq. (|TT1) and expresse 



F(y,z) = / j^^e^'^y-^^FLii) + Fs 
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each 1/ into 1/l + 1/s to obtain 



Since (f^)^ = = (fs)^; the above equation then becomes 



P P' 



P 



P 



(12) 



(13) 



where the first term fifs in the right hand side leads to a vanishing result as it contacts with 
the term i^/P term. The only contribution can only come from the second term in the 
right hand side of Eq. (JT3l) . In addition, the ^ can be expanded in the terms proportional 
to ^, fi, This gives 



In ■ I 



Due to = = 0, it further reduces to 



In- 1 



By substituting the above back into Eq. ( 1121) . Eq. ( jlTj) is then obtained by noting that 

Using Eq. ( ITTl) . one can systematically include the effects from the non-coUinearity and the 
off-shellness of the coUinear partons. This property of the long distance part of the parton 
propagator plays an important role in our following analysis, and its meaning will become 
more clear after we have investigated real cases lately. 

According to the parton model, the hadron amplitudes are defined as the probability for 
finding the on-shell partons inside the hadron. The parton amplitudes are then required to 
contain only the on-shell components of the external parton momenta. However, according to 
Eq. ([8]), either the on-shell momentum or the collinear momentum / can be assigned for an 



on-shell parton. 



on-sn eii p art 

0, y, H [n 



factorization 
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'herefore, there arise two factorization schemes, the collinear factorization 

and t^ 



ISjfQCD factorization 
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the /ct factorization Q, H, Q, Q (PQCD 
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28(1). In the fc-r factorization scheme, an on- 



shell parton carries a momentum Z^. On the other hand, in the collinear factorization scheme, 
an on-shell parton carries a momentum /. In this work, we follow the QCD factorization 
to use the collinear factorization scheme as our basics. Our proposed collinear expansion 
method is composed of following steps: 



1. Use scale analysis for the parton amplitudes according to the scales of parton momenta 
to find out the leading regions of the parton momentum configuration. 

2. The parton amplitudes are expanded into a Taylor series with respect to the leading 
regions of parton momenta. 

3. The expanded parton amplitudes are substituted back into the contraction with the 
hadron amplitudes to extract relevant contributions up to specific twist order. 

4. The factorization of parton momentum integrals is accomplished by means of integral 
transformations (See, for example, Eq. fl2^ ). 

5. The color structure of the parton amplitude is extracted to be attributed to the hadron 
amplitudes to complete the color factorization. 

6. The factorization of spin indices is completed by means of Fierz transformation. 

7. The property of the long distance parton propagator is used to extract higher twist 
contributions. 

We are now ready to discuss the collinear expansion. First, we order the parton ampli- 
tudes in as 

Tl,{l)=Tl^+Tl^{J) + 0{a^, (14) 
^1.(^1, ^2) = Tl^di, h) + Tl^d,, h) + 0{al) , (15) 
Ti\lB, ^Mi5 h'h) = tI^'^''\Ib, Imi, h'h) + 0{al) i (16) 

where the superscription (0) and (1) are used to denote the relevant parton amplitude of 
zeroth and first order in a^, respectively. The parton amplitude T^j^^ is just the tree vertex 
TiSij in the diagram as depicted in Fig. [21 There are vertex and penguin diagrams for T^^*^^^ 
amplitudes as depicted in Fig. HI The parton amplitude T^j'^J describes the tree diagrams 
as depicted in Fig. [3], in which two quark partons and one gluon parton from the meson M2 
are interacting with a local four fermion operator Oi = (^iFjfe) (52^1^3) with rj(ri) the Dirac 
gamma matrix for the operator. The parton amplitude T/^ starts from 0{as) diagrams as 
depicted in Fig. O and is denoted as T^^^^\ 
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We first expand the parton amplitudes T^/^ with respect to the colhnear components of 



their relevant parton momenta as 




2 r^rpim 



where 




I 



(18) 



(19) 



The expansion series are then substituted back into the convolution integrals with the hadron 
amplitudes for further analysis. The reason why one can expand the parton amplitudes with 
respect to the relevant collinear momenta will be explained in detail below. 

A. Expansion with T-^^°^ 

The expression for the contributions associated with T^^^'^ is written as 



where the loop momentum I is carried by the loop parton in Fig. [T](a). Since the parton 
amplitude T-^*'"'* is independent of /, we propose to use the following integral identity to 
transform the expression into a form consistent with the parton model picture 




(20) 




1 . 



(21) 



The transformed result appears as 




(22) 



where 




(23) 
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The following integral transformation has been used in the above to simplify the expression 

dS5^'\y-\n)G{y,Q) 



= G(An, 0) , (24) 

where G{y, 0) denotes any function of the coordinates. Two comments for the above integral 
transformations Eqs. (12T|) and (l24l) are necessary. First, the momentum fraction x for the 
parton of the meson M2 is introduced. Second, the quark field q{Xn) is ordered in light-cone 
direction n. This implies that the hadron amplitude ^m2{x) is defined on the light-cone 
= where n'^ is a null light-cone vector. By using the above integral transformations, the 
parton amplitude and the hadron amplitude are only related by the momentum fraction x. 
Because the parton amplitude T-^''°^ is equal to Fj^jj, the integral over x is then associated 
with $A/2(x). 

The factorization of spin indices depends on the structure of F,. For {V — A){V ± A) 
operators, Ti = (V ± A) and can be expanded into ^(1 ± 75), ± 75) and 7^75. If M2 is 
a pseudo-scalar meson, only the axial vector part can contribute. However, only 1^75 leads 
to leading twist contributions. The ^75 will result in twist-4 contributions and 7_l75 will 
not contribute. For other types of meson, similar considerations can be made. We now 
explain how the ^75 part can contribute. The long distance part of the parton propagator 
can interact with vertex ^75 to have 

It is also applicable for the other parton propagator of the anti-quark line. The short distance 
part of the parton propagator ifi/{2n ■ I) and the vertex ija are absorbed into the parton 
amplitude. This results in 

»i 

dxTY[Tij^ap{x, X, x)w'^,Wp,^lj'^Q{x, X, x)] (26) 





where w^, = g'^, — q^Ua' , 

2n-l 2n-l 



Tij^ap{x,x,x) = (i7a)7r^^?°^7r-^(-^7/3) 



<f;3(x,a:,x) = ^e^^^(M2|g(An)^9"'(An)z9^'(An)g(0)|0). (27) 
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There are corresponding contributions from the two gluon insertion diagrams depicted in 
Fig. [6|, whose expression is written as 

/ (^I'^'I^S ('.. h)n,>l,Kjlu fe. '3)1 (28) 

where we have employed the hght-cone gauge n ■ A = for the gluon fields, and the parton 
amplitude and hadron amplitude are expressed as 



x(M2|g(|/)(-^7A"(z))(-(7A^H)g(0)|0) (29) 

Since T-^^^p{h,l2,h) can be replaced by T^j^p{xi,X2,X3) straightforwardly, the momentum 
integrations over li,l2,h can be transformed into the integrations over xi,X2,X3. We then 
obtain 

j dXidX2dx3TT[T.I^Jf^{xi,X2,X3)^"j^^^^{xi,X2,X3)] . (30) 

The combination of Eq. (!26l) and Eq. (!30|) gives 

j dxidx2dx3Tr[T.j"Jp{xi,X2,X3)w2,w'^,^l/^ j^{xi,X2,X3)] (31) 

where 

^ti^{x,,X2,X3) = J^^J^J ^e^^-e-(---)e-(---) 

x{M2\q{Xn){tD''{7]n)){iD^{ujn))q{0)\0) (32) 

with iD'^ = id°'—gA°' being the covariant derivative. Since fi is of 0{Q~^), T^j^^ is of 0{Q~^) 
as the scale of T^j^^^ being of 0(1). The relevant contributions are of higher than twist-4. 
The above example is to show that, using the coUinear expansion, one can calculate the tree 
level higher twist corrections from the dynamical partons in a systematic way. Because we 
only intend to calculate the twist-3 corrections, we will not further explore the contributions 
of twist order higher than three. For —2(5 — P){S + P) operators, fj = 75. Up to twist-3, 
the expression appears as 

I\xTt[T;I'^<!>m,{x)] = -l^'cia;Tr[7;fV75]IV[ri75$M,(x)] 

+ ^ rfxTr[7;f )75]Tr[75$A./,(x)] . (33) 
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By identifying Trlfi'y^^ Mii^)] and Tr[75$M2(3^)] as the twist-2 and twist-3 two parton LCDAs 
of the M2 meson 

Tr['i/ij,<^>M,{x)] = -^/M.0it(^) ' (34) 
1V[75^M.(x)] = -ifM,f^x<pMlp{x) , (35) 

where fi^^ = m\[^ / {fhq + ffiq) with fhq and fhq the current quark masses and mu^ the meson 
mass, we recover the naive factorization result up to twist-3 order. 



B. Expansion with T^^^^ 

''Jit*' 

We show the hght-cone gauge n ■ A — and the covariant gauge 9 • ^ = for the 
expansion with T^^^^ . Because the analysis is tedious, we outline the procedure, here, and 
leave the details for next section. The first step is to take a power counting for the parton 
amplitude T^^'^^ . There are three interesting regions. The first region is composed of either 
two soft loop parton momenta, or one soft loop parton momentum and one coUinear loop 



parton momentum. The T^^^^ in the first region is counted as A~^. The second region i 
composed of two coUinear loop parton momenta. The T^^^^ in the second region is counted 
as QX"^. The third region is composed of either one coUinear loop parton momentum and 
one hard loop parton momentum, or two hard loop parton momenta. In the third region, 
the tI^^^ is counted as Q~^. We conclude that the region composed of two coUinear loop 
parton momenta is dominant. 

Let's, first, consider the light-cone gauge n ■ A = Q. The expansion series of T^^^^ili^ I2) 
with respect to Zj, i = 1, 2, are written as 

Tl^Hu h) = Tl^iii, k) + ^^^^ ^"2)(/. - IkY + • • • , (36) 

k=l,2 

where Tl^^^^^{ii, 4, ^2) are defined by assuming the low energy theories 



IS 



(37) 

The expansion series are then substituted back into the convolution integrals. Since the 
gauge condition n • A — with a light-cone vector satisfying n • q — 1, — 0, and 
n- 1± — 0, the first term ^^^["'(/i, ^2) leads to the result 

/ / (£^T^[^- ^^^)-«2(^i, h)] (38) 
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where we have introduced = Oa ~ q^Ua- Similarly, we employ the following integral 
identities to simplify the momentum integrations 

Jo Jo Jo 27r Jo 271 
The expression appears as 

/ dx, I dx2Tr[7;^(;)(a;i,X2X$^,,(xi,a;2)] (40) 
Jo Jo 



in which 

= T^j^iliM ^ , , (41) 

Jo Jo 2TT 

x{M2\q{Xn){-gA^{7]n))q{0)\0) . (42) 

In the above equations, we have used the following transformation for any function G{y, z, 0) 
of coordinates y'^ and 

d% r d% 

= G{Xn,r]n,0) . (43) 

For tree diagrams, only color singlet operators can contribute at twist-3. The remaining 
task is to finish the factorization of spin indices. For {V — A){V ±A) operators , the related 
contributions are of twist-4, which are beyond our accuracy. For —2{S—P){S+P) operators, 
the related contributions are of twist-3. To obtain the collinear limit k Xi, i = 1,2, for 
tI^^^{xi,X2) , we have used the following substitution for the parton propagators 

i - ^ . (44) 

This is because the vertex 27^ in T^^^^J{xi,X2) is transversal and only the off-shell part of 
the parton propagators can contribute. See further explanations in the next section. The 
terms associated with T^jf!^^^, are of twist-4 and higher. They are neglected accordingly. The 
factorization of spin indices results in 

I fo r ^^2Tr[7;^2V„;375]<Tr[a"^5$^t;,(xi, X2)] . (45) 
The explicit expression for Eq. (j45|l is left to the next section. 
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We now consider the expansion with covariant gauge d ■ A = 0. Since the factorizations 
of the momentum integrals and color indices are independent of gauge condition, we can 
go through to consider the factorization of spin indices. The first term in the expansion 



of T^jJ{xi,X2), under covariant gauge, is related to the gauge invariant phase factor of 



the related two parton amplitudes. The gluon fields A^^ in ^^^^ can be expanded as A^^ 
n ■ Aq^ + q ■ An'^ + d'^A°'. The contraction TtItII'^J {xi, X2)q'^n ■ $^2(2^1! 2:2)] leads to 



J dxi J dx2TT[Ti^^^-^n ■ ^M2{xi,X2)] , (46) 
where k = I2 — h being the gluon momentum and 

n-^M,{xi,X2) = / — / M'^^e'^^^^-^^HM2\qiXn){-gn-A{rjnM0)\0) . (47) 
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The terms with q ■ An'^ vanish since the covariant gauge condition d- A = 0. The terms with 



the contraction TilT^j'^J {xi, X2)d'^^1.j^{xi, X2)] are of higher twist than twist-3 and will be 



neglected. With the above considerations, the contraction Tr[T^^.^°''(xi, X2)$a^2(^1' -^2)] leads 
to contributions of twist-2 or higher than twist-3. 



We next consider the contraction Ti\rl^^\^{xi,X2){lk — 4)'^'^'mo(^1' ^2)], which can be 



rewritten as 

Tr[^S,.(a;i, Xfe, X2X,$aJ;(xi, x^, X2)] (48) 

with 

$£(xi,Xfc,X2) = /^/ ^e^'^^e^'^(^-^^)(M2|g2(An)^^7G'^'^(r/n)g3(0)|0) . (49) 

Note that only transversal part ^(4 — Ik)^ of the {Ik — IkY can contribute at twist-3. For 
(y — A)(y ±A) operators, the contributions are of twist-4. For —2{S — P){S + P) operators, 
the result appears as 

x{6{xk - xi) + 6{xk - X2)) . (50) 

The reader may have noticed that the terms in the expansion series of T^^^^J{li, I2) in Eq. fl36|) 
are of different twist order under the covariant or the light-cone gauge. For example, the 
twist-3 contributions are from the the first term in the expansion series of T^j'^^{li, I2) under 
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the light-cone gauge. On the other hand, under the covariant gauge, the twist-3 contribu- 
tions are from the the second term in the expansion series of T^^^^J {li,l2). Since the parton 
amphtude and hadron amphtude under the colhnear expansion are required to be gauge in- 
variant, respectively, this feature of the collinear expansion method can be used as a guiding 
principle for calculations. 



III. TWIST-3 CORRECTIONS 



In this section, we make a more detail descriptions for the twist-3 contributions from the 
three parton Fock state qqg of the M2 meson. The amplitude for the three parton qqg of 
M2 interacting with the operator Oq at the tree level for 13 — > M1M2 decays is written as 

(Mi|gi(0)(l-75)K0)|5) 

+ (7-f+f)t+,, (+'gA(y))(i + 75)]g3(0)|0) . (51) 

The / and k denote the momenta carried by the q2 quark and g gluon fields in Fig. [3](a) and 
(b). We first employ the light-cone gauge n ■ A{y) = 0. The gluonic fields A°'{y) represents 
A°''°' {y)T"' with the color matrix in the fundamental representation '^t°'T'' = 5""^ /2. To 
relate to the previous introduced collinear expansion, we recast the convolution integration 
part of Eq. (15T1) into the form 

/ (|7(0iTr[T„'<»)(fc,i)<.*'"(*.Ol (52) 

where the parton amplitude T/i^^\k,l) is defined as 

^.""'(M) ^ lin.)^jf0- + ^:}t;l;^A -nM^ + 7.) (53) 

and the meson amplitude ^^'(k,l) 

^^'{k,l) = j d'y j d^ze''-^e'^<y-^\M2\q2{z){-gA>^\y))q,{S)m . (54) 

The tensor w^, = g^, — q^n^, has been introduced. Note that, for convenience, we have made 
a change of variables for the loop parton momenta I = and k^ = {I2 — li)^ ■ We assume 
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that the emitted M2 meson is highly energetic. As shown in last section, the dominante con- 
figuration is composed of coUinear li and I2. This allows us to expand the parton amplitude 
"with respect to / = xq and k = {x' — x)q 

(/ - ly + '''' 



l=Lk=k 



{k-ky + ■■■ (.55) 

i=ih=k 



Substituting the first term back into the convolution integrations gives 

j dx j dx'i:Y[Tl^^{x,x')w^^,^^' {x,x')] (56) 



where 



d^l f d^k 

<l>^'(x',a;)= / 6{x - I -n) / -——-5{x'-x-k-n) I d'^y / d^z 



(2Txy ' ' J (27r) 

yie^'-^e^'y^M^Mzy-gA^'iyyq.mQ) ■ (57) 

In the above coUinear limit step T^j^J{l, k) T^^^^J{x, x'), there arises an infrared divergence 
as x' — > 0, which is from the denominators of virtual quark propagators 

(x'qy + te ■ ^ ^ 

We regularize this divergence by the following method. Since the full quark propagator with 
momentum V = I + k can be decomposed into its long distance part and short distance part 
as 



{l + kf + ie (I'Y + ie 2n-l' 

The long distance part gives vanishing result upto twist-3. The short distance part is 
absorbed by the parton amplitude. The divergence is then re gula rized by replacing the 
quark propagators with its corresponding special propagators [zl. \a\ 

ix'^ ifi x' — X 

[x'qy + ie 2x + ie x' — X + ie 

The introduction of a special propagator for an on-shell fermion propagator is due to the 
fact that the fermion propagators in Fig. 3(a) and 3(b) become on-shell and divergent after 
the coUinear expansion. The divergent part of these propagators leads to long distance 
contributions that should be included into the twist-2 distribution amplitude for the M2 
meson. However, there are also finite contact part of these propagators, which leads to 
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contributions of one twist higher. The more detailed explanation about the meaning of the 



special propagator refers to 



7 



4Q| 



Under light-cone gauge n ■ A{y) = 0, it is convenient to transform the gluon fields A^{y) 
into its field strength G^^{y) by using following replacement 

and 

<(^>^')-^<(^,^')- (62) 
The factor iny/{x' — x) is then absorbed by tI}'^\x,x') into the form 



T,?)(x,x')=T;(°)(x,x')^. (63) 



The factorization of the spin indices gives 

IJdxJ cix'Tr[T^^W(a;,x')a„^75]<TrK^5$^^'(x,x')] + ■ ■ ■ , (64) 

in which other spin decompositions give higher twist contributions. 

The numerators in the contraction Tt[T^I)^^ aai3'y5\ can give terms proportional to 
nun^iqarip-riaqfs), n^d±^^^"{qanf} - riaq/s) , and n^d_i_^^^»e±^ai3- The transversal tensors rf_L,a/3 
and e±,a/3 are defined as d±^af3 = qa^p + qpria — Qap and e^^^ap = ^apr^xq^n^- The trace of 
dj^^ap is defined to be negative d°[^ = —2. Since v and yU indices in are antisymmetric 
under /i z/, the terms proportional to ni,n^{qanf3 — Uaqp) then vanish. For those terms 
proportional to n,yd±^^^ii{qanfs — Uaqp), as they are contracted with Tr[cr"^75$''^'(x, a;')], the 
qa factor in riyd^^^^ni^q^np — Uaqp) results in twist-4 contributions by using the property of 
the long distance propagator of the quark fields. The terms proportional to n^d±^^^iie±^ag 
lead to twist-3 contributions. The final result appears as 

f G'^(x,x')nfn0 ^ ^ 

dx / dx' 7 \ ^ , 65 
J [x' — x)x 

where the function C^^x, x') is defined as 

= j -^J^^i^-^-^) J -^^i^i^' - ^ - ^ ■ ^) j A j d^z 

xe''-^e'^-y{M2\q2{z)a,^^^wl,gG^'''{y)qM\Q) . (66) 
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Note that we have used the G-parity symmetry x ^ x' to simphfy the above result. This 
assumption is valid for vr mesons, but may not be appropriate for the K or rj mesons. 
Therefore, it is noted that, in the above result, there exist symmetry breaking effects for K 
and 1] mesons. However, we will ignore such a corrections from the symmetry breaking in 



the following calculations. By referring to the definition 



3 



{M2\q2{z)(r^,-f5gGa/3{y)q3mO) 

iqaqi,d±,up - qaQyd^^^p - qf3q^d±^^a + qf3qud±,ap,)Tiz, H , (67) 



where 



T{z,y)= dx (ix V^-"e-'("-"')''-^T(x, x') , (68) 
Jo Jo 



we can arrive at the result 



J {x' — x)x 

^■'h.,m, [ ,^.n^ (69) 



By using the normalization for (Mi|gi(0)(l — 75)6(0) it is easy to derive the tree level 
three parton contributions for operator Og as 

{0,),^,luon= '^"'7"- , (0l)/ (70) 

with 

(71) 

We now explain the expansion with the covariant gauge d ■ A = Q. We first decompose 
A'^[y) into its longitudinal and transversal components as A^{y) = n ■ A(y)q^ + d'^ ^,A^^' (y). 
The transversal part d'^^,A'^'(y) results in contributions of higher than twist-3 . The lon- 
gitudinal part n ■ A{y)q°' gives twist-3 contributions. Similar to the light-cone gauge, we 
need to transform the gluon fields into its field strength. Here, it needs one transversal 
momentum k± factor from expansion of the parton amplitude Tl}^\l,k) in Eq. ( !55l) . The 
contraction of tI}^\x, x') with q^n ■ ^m2{x, x') leads to two parton gauge phase factor 

7-/(0)/ A _ 7-/(0) ('^^ 

Tr[r;(°)(x,x')g'^n ■ <fM.(x,x')] = Tr[ ■ $M.(a;,x')] (72) 
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It is convenient to write {k — ky = d^j_p,{k — ky' + q ■ kn^. Only transversal part k'^ = 
d±p'{k — k)p' contributes at twist-3. This is because the term OT/i^'^^ / dk'^ can have terms 
proportional to Qf^u and cx^u- The terms related to q ■ kn" leads to twist-4 contributions. 

For the transversal part k±, only a^^^, terms can contribute. Let the fc^ factor absorbed 
into $'^(/, k) and use the replacement k'^A^{y) —iG'^^{y), we can derive the result 

x(Mi|gi(0)(l-75)K0)|5) , (73) 



where 



and 



dk- {x'-x)xq^^ 



G'^ix, x') = I 7^5(x -l-n) [ 7^5(x' - x - k ■ n) I d^y I d^z 



(27r)4 ^ 'J (27r)4' 

xe''-^e'^-y{M2\q2{z)igG^^'{y)q,Xm) ■ (75) 

The contraction of cr^,y with G"^^(x, x') gives 

^r[^a,^G^>^{x, x')] = Z^!^^!!^t(x, x') . (76) 

Note that the q^ factor in the denominator of dT^^^\x' , x) / dk^ is cancelled by the q^ factor 
in the numerator of Tr [zcTj^^G"^'^ (x, x')]. It is easy to see that Eq. fl73|) is equal to the result 
derived from the light-cone gauge. This explicitly shows the gauge invariance of the three 
parton contributions. 

There are related diagrams, such as those in Fig. [3](c) and[21^d). Because the spectator 
quark of the B meson can carry only a soft momentum, this makes the relevant contributions 
associated with Fig. [3]^c) and [3](d) dominated by soft gluons as the form factors . 
In addition, the relevant contributions are of 0(777^^) with respect to the leading twist 
amplitude. It can be understood as following. The sum of the lower parts of the diagarms 
in Fig. [3](c) and [3](d) is proportional to 



(Milgi 



2p^ + -f4 _ 2Pf,^ 

-L 7 i 7 



2p-k 2Pb-k 



b\B) , (77) 
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where k is the momentum of the gluon from M2 and the equation of motions for b and qi 
quarks have been used. After taking the coUinear hmit, k — > x'q, we write the expression as 

Aq, + B^,q^ , (78) 

where 

A ^ -l^{M^\q^T,\B) , (79) 



x'Q2 
1 



= 7^77^{Mi\qih.l^T, + ra>.7.)b\B) . (80) 



For hght-cone gauge, only A term contributes. As for the covariant gauge, only B^^, term 
contributes. The only contributions come from {V — A)(y±A) operators. This implies that 
the upper parts of the diagarms in Fig. [3]^c) and[3]^d) are proportional to the twist-4 LCDA 
of M2. The combination of the upper and the lower parts gives a 0(m^^) contributions with 
respect to the leading twist amplitude. 

There are possibilities that the additional gluon of the M2 meson can interact with the 
spectator quark of the B meson. Since the spectator quark carries a soft momentum, the 
momentum conservation at the interaction vertex prevents the momentum of the gluon from 
being coUinear to the M2 meson's momentum. Therefore, there require additional radiative 
gluons interacting between the other parton lines and the spectator quark line to make the 
momentum of the gluon to be coUinear to the M2 meson's momentum. This results in 
contributions of order 0{as)- We identify the relevant contributions as 0{as) three parton 
corrections. As mentioned previously, we plan to discuss these contributions in other places 



44|. 



The total twist-3 contribution from operator Oq is then equal to 



For operator Og, there are similar results. 
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IV. APPLICATIONS 



For penguin dominant B ttK decays, the relevant decay amplitudes under QCD 
factorization are parametrized as the following 



(«4 - \aio)+r^K - ^«8) 



A 



ttK 



-V2A{B- ^ n'K-) = [A„ai + Ap« + + V?(«6 + O]^ 

3 

+ [Xua2 + Ap-(-a7 + ac))]AKn 

+ i>^ub2 + {\uXc){h + bf'^))B^K , 

-A{B^ ^ n+K-) = [A,ai + Ap« + <o) + + alM^K 

1. 



wK 



+((A„ + A,)(63- 

V2A(B^ TiK^) = A{B- n-R^) + V2A{B- n^K') - A{B^ tt+R-) (82) 

where Ap = V^bV^*.., ca = aii^iiK), and Apof = Xua^ + AcO^. The CP conjugation of decay 
amplitudes are obtained by replacing Ap A* for the above amplitudes. The factorized 
matrix elements are defined as 

Ak. = ^^{ml - ml)Fi-^{ml)U • (83) 



The form factors are defined 



2 2 

m R — m 



{P{p)\qrb\B) = F^-^{q'){P^ + + [Fo^^^(g') - F^^^ {q')] "'^ q'^ . (84) 

The form factors coincide as = 0, F^-*^{0) = F^^^^^O). The expressions for the param- 
eters are referred to For numerical calculations, we will use the following input 
parameters 

A^^g = 0.225GeV , mfe(mb) = 4.2GeV m^{mb) = 1.3GeV , m,(2GeV) = 0.090GeV , 

iKfel = 0.41 , \Vub/V,b\ = 0.09 , 7 = 70° , t{B-) = 1.67(ps) , ^^^^ 

(85) 

T{Bd) = 1.54(ps) , U = 131MeV , fx = 160MeV , fe = 200MeV , 

F(f-^ = 0.28 , Fq^^^' = 0.34 . 
For A„ and Ac, we take the following convention for their parametrization 

^ = tan^ e^be-''^ (86) 
Ac 
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where 

tan^ ^ 
Rh 



1- A2 ' 
1-AV2 Kb, 
A WJ 



A = IK.I . (87) 

The value of A is taken as 0.22. 

By using previous input parameters, we hst the values of Oj, i = 1, • • • , 10, and bj, 
j = 1, ■ ■ ■ ,3, calculated at the scale rrih = 4.2GeV as below 

ai = 0.995 + 0.018Z , 03 = 0.209 - 0.104z , 03 = -0.003 + O.OOSi , 

= -0.031 - O.OlSi = -0.030 + 0.027i , 05 = 0.007 - 0.004i , 



r^a^ = -0.050 - O.OlSz , r^ag = -0.047 - O.OOSz , ay/a = 0.007 + 0.006z , 
r^a^/a = 0.087 - 0.043i , r^al/a = 0.094 - 0.021i , ag/a = -1.135 - 0.024i 
a«o/a = -0.175 + 0.093z , alja = -0.175 + 0.093i , r^6i = 0.021 , 
r^62 = -0.008 , r^fes = -0.006 , r^&f ^/a = -0.018 , 



(88) 



in which 



" A^K " m|Fo^-(0) ' ^^^^ 



and 



2m 



2 



= f 1 , ■ (90) 

mb{mq + ms) 

The tree level three parton contributions modify the parameters as r^(l + A^f^) with 
parameter A^^ defined in Eq. f l7ip . The value of A^_^ depends on the model of the three 
parton distribution amplitude Tm2{x, x'). Here we employ the model derived from the light- 



cone sum rule 



4l| 



T{x, x') = 360rixx'{x - x'f{l + ^{7{x - x') - 3)), (91) 

where the parameters are assumed to be 77 = 0.015 and u = —3.0 for M2 = vr, K, or 77. This 
give us A^l = 0.585. 

The branching ratio for a, B ^ ttK decay is given by this expressions 

Br(B^7iK)= \A(B^7iK)\^. (92) 

167rms 
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We can use the above formula to predict CP averaged branching ratios for B irK decays. 
The predictions with three parton corrections in units of 10~^ are given as 



Br{B- - 




= 19.0 , 




Br{B- - 




= 10.0 , 




Br{B^ - 


TT+K-] 


= 16.1 , 




Br{B^ - 




= 7.7 . 


(93) 



For comparison, we also list the predictions with only two parton contributions in units of 
10^^ in the following, 



Br{B- - 




= 11.2 , 




Br{B- - 


n'K-] 


= 6.1 , 




Br{B^ - 


n+K-] 


= 9.4 , 




Br{B^ - 




= 4.4 . 


(94) 



For reference, we enlist the experimental data in units of 10 ^ summarized by the HFAG 



group 
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Br{B- n~K^) = 23.1 ± 1.0 , 
Br{B- n^K^) = 12.8 ± 0.6 
Br{B^ n+K^) = 19.7 ± 0.6 

5r(fiO ^ 7r°i^o) = 10.0 ± 0.6 . (95) 

By comparing the predictions with or without the three parton corrections, one may notice 
that the the predicted branching ratios are significantly enhanced by about 1.65 ~ 1.75 
times in their magnitudes. 

The two parton predictions for the B ttK decays made here are consistent with the 
findings of previous studies using QCD factorization approach js], I43I. The two parton 
predictions made in |43] are much lower than the experimental data under the QCD factor- 
ization approach. Because the calculations of the three parton corrections were inaccessible 
in their studies, this led them to conclude that the QCD factorization is impossible to ex- 
plain the penguin dominant B ttK decays. The two parton predictions made in |^ for 
B ttK decays are still lower than the data. Only extending the predictions by using 
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extreme limits of input parameters can make the predictions to be consistent with the mea- 
surements. This seems not a reasonable solution from the theoretical point of view. On the 
other hand, as shown in the above, our approach has shown that the predictions with three 
parton contributions are more close to the data than the two parton predictions. Although 
the predictions with three parton contributions are still lower than the experimental data, 
the 0{as) corrections can improve the predictions. 



V. DISCUSSIONS AND CONCLUSIONS 



The significance of the three parton contributions for penguin dominant B ttK decays 
can also be seen from a phenomenological point of view. By appropriate arrangement, the 
parameters can be calculated for B -kit decays. For the pure penguin B^ —>■ n^K^, 
its dominant contributions arise from the \aK7cK) + r^ag(7ri^')| term. The uncertainty due 
to the form factors can be eliminated by considering the ratio between the decay rates of 
B^ Ti^K^ and B^ tt^tt^ as the following 



aifvTTr) + a2[TTTc) 



\Vub\ 


fn 


'T{B- -R-K^) 


\Vcb\ 


fK 


_2T{B- TT-vrO) 



1/2 



0.105 ±0.001 



(96) 



where the error comes from the branching ratios. In the above, we have used the branching 



ratio BriB 



TT TT 



5.7 ± 0.4. According to QCD factorization calculations, |ai(7r7r) + 



a2(7r7r)| = 1.17 and \a^^{T{K) + r^ag(7rir)| = 0.08 with only two parton contributions. This 
gives the prediction of the ratio to be 0.066, which is lower than the experimental value. By 
adding the tree level three parton contributions, the factor then becomes ^"^(1 + ^/^'^) 
and makes \a^^{TiK) + r^i)- + A^)a'^{'nK)\ = 0.104. The predicted ratio is changed to be 
0.089, which is closer to the two parton prediction. The above fact may also indicate that 
the three parton corrections could be important for understanding the penguin dominant 
B ttK decays under the QCD factorization approach. 

In order to make sure that the three parton contributions are indeed significant and also 
compatible with the QCD factorization, an important task is to finish 0{as) calculations 
for the three parton contributions. The related work in this direction has been proceeded 



and will be reported in our another preparing paper 



4^. 
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There are similar three parton contributions having been calculated under the light- 



cone sum rule 



45 



461]. In 45(1, the twist-3 three parton contributions associated with soft 



gluons are calculated in the framework of light-cone sum rule. The contributions are shown 



negligible in B —>■ tttt decays. In 



461] , the contributions from the diagrams similar to Fig. [31(c) 



and [21(d) were calculated for B ttuj decays. They are found to be twist-4 and vanishing. 
In addition, significant effects were found due to the three parton Fock state of the vr in the 
B -Kuj decays. Since they are dominated by soft gluons, it is better determined by QCD 
sum rule. As shown in {46], in the Euclidean region of {p + g)^, the relevant contributions 
are from the twist-3 and twist-4 three parton LCDAs of the vr. As mentioned before, we 
identify these power corrections as non-partonic ones. From the theoretical point of view, 
we suggest that the partonic and non-partonic power corrections should be distinquished 
under the QCD factorization, although they may be equally important in phenomenology. 

For comparison, we employ the replacing rules for the Oj coefficients 46|] to account for 
the three parton effects from the Mi meson. The rule is 



a2i ^ a^i + [1 + {-lY''^''^C2i-ih/2 



(97) 



where i = 1, ■ ■ ■ ,5, and Ci are the Wilson coefficients calculated at the scale ^ih = 1-45 GeV, 
and /s = 0.12, which is assumed to be universal. With these three parton corrections, the 
predicted CP averaged branching ratios for B ttK in units of 10~^ are 



Br{B- 
Br{B- 
Br{B^ 
Br{B^ 



9.5 



. n^R-) = 5.9 , 
n+R-) = 8.4 , 
n^RO) = 3.4 , 



(98) 



which becomes smaller than those predictions in Eq. (p^ . 
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FIG. 2: The Feynman diagram describes the tree level four parton amplitude, T^j . The square 
symbol represents the vertex of weak interactions. 







(a) 



(b) 



(c) 



(d) 



FIG. 3: The Feynman diagrams describe the tree level five parton amplitude, T^j^'^J . The square 
symbol represents the vertex of weak interactions. 
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(a) (b) (c) (d) 




FIG. 4: The Feynman diagrams describe the 0{as) four parton amplitude, T-j ' . The square 
symbol represents the vertex of weak interactions. 




(e) (f) 



FIG. 5: The Feynman diagrams for the 0{as) six parton amplitude, T^^^^\ The square symbol 
represents the vertex of weak interactions. 
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FIG. 6: The Feynman diagram for the six parton ampUtude with \qqgg) Fock state. The square 
symbol represents the vertex of weak interactions. 
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